Recent developments in POWHEG 



Paolo Nason* 

INFN, Sez. di Milano Bicocca 



E-mail: Paolo . Nason@mib . inf n . it 



I review recent developments in POWHEG, a method for interfacing parton-shower generators 
with NLO QCD computations. I illustrate recent progress in understanding several features of 
the method, and in clarifying similarity and differences with respect to MC@NLO. Furthermore, I 
briefly describe a recently introduced framework, the POWHEG BOX, that allows the automatic 
POWHEG implementation of any given NLO calculation, and has been recently applied to Z+jet 
production and to Higgs production via vector-boson fusion. 



RADCOR 2009 - 9th International Symposium on Radiative Corrections (Applications of Quantum Field 
Theory to Phenomenology) , 
October 25 -30 2009 
Ascona, Switzerland 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Recent developments in POWHEG 



Paolo Nason 



The POWHEG method is a prescription for interfacing NLO calculations with parton shower 
generators. It was first suggested in ref. [|]], and was described in great detail in ref. [^. Until now, 
the POWHEG method has been applied to Z pair hadroproduction [||], heavy-flavour production [Q], 
e~^e^ annihilation into hadrons [^] and into top pairs Drell-Yan vector boson production |0|8|], 
W' production [^], Higgs boson production via gluon fusion [|l^, |ll|], Higgs boson production 
associated with a vector boson (Higgs-strahlung) |[IT|], single-top production [ p^ ] Z + 1 jet pro- 
duction [13], and, very recently, Higgs production in vector boson fusion [^]. Unlike MC@NLO 
[p4[|, POWHEG produces events with positive (constant) weight, and, furthermore, does not de- 
pend on the subsequent shower Monte Carlo program. It can be easily interfaced to any modem 
shower generators and, in fact, it has been interfaced to HERWIG [15, 16] and PYTHIA [17] in 
refs. [||0,[Tg,|T2|]. 



1. The concept 

The basic concept of POWHEG is better clarified by considering a process with a single mass- 
less coloured parton involved (one can consider, for example, semileptonic top decay, neglecting 
the b mass). In a shower Monte Carlo, the radiation of a final state light parton is generated with 
an algorithm that resums all leading log corrections to the Born process. The hardest emission in a 
shower Monte Carlo is generated according to the formula [^ 



da = B(<t>B)d<t>B 
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where t is the radiation transverse momentum, to is the minimum allowed value for t (typically of 
the order of a hadronic scale), B{(^B)d^B is the Born differential cross section, and E^^d^^B^r is 
the real radiation differential cross section in the Monte Carlo (MC from now on) approximation. 
It is assumed that the full phase space including radiation is parametrized in term of the Born phase 
space and the radiation phase space <I>r, i.e. = 0(<I>b,<I>^). In typical MC's, the radiation 
phase space is determined by three variables characterizing the collinear splitting process, like, for 
example, the splitting angle, the momentum fraction and the azimuth. The radiation transverse 
momentum ? is a function of and O,.. It can be defined as the momentum component of the 
radiated parton orthogonal to the momentum of the radiating parton. The MC Sudakov form factor 
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represents the probability for not having radiation harder than 

The basic idea in POWHEG is to improve formula < \l.l\ ) in such a way that NLO accuracy is 



reached. One replaces formula (1.1) with the following one 
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where is the full phase space, with d<t> = d^Bd^r, and R is the exact radiation cross section. We 
require that /? — )■ 7?^ in the soft and collinear limit, and that R" < R, so that the last contribution in 
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the square bracket of eq. (1.3) is non-negative. The choice R^ = Ris also possible, and it is quite 
common. We have also defined 



V(<^b)+ R'md<^r 



(1.4) 



where V is the virtual NLO correction to the Bom process. Notice that soft and coUinear singular- 
ities in V cancel against those arising from the integral of R^ in the square bracket of eq. (|]^. The 
Sudakov form factor is now 



exp 
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The phase space factorization in the POWHEG formula needs not to match that of any shower Monte 
Carlo. One only requires that in the soft and collinear limit the full phase space is related to the 
Born phase space <I>5 in the correct way, i.e. they are identical in the soft limit once the soft particle 
is removed, and they are identical in the collinear limit once the collinear particles are merged. The 
POWHEG formula ( [L3| ) can be viewed as an improvement of the Monte Carlo formula ([Ell), such 
that the Bom cross section is replaced with an NLO inclusive cross section, and high transverse 
momentum radiation is corrected so that it becomes exact at large angles. In fact, for large t the 
Sudakov form factor becomes 1, and the POWHEG cross section reduces to 



do = Bx —d^+\R 
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(1.6) 



since B/B = 1 + ^{(Xs). At small t the POWHEG formula becomes equal to that of a standard shower 
MC, up to higher order terms. However, since by construction 
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the POWHEG formula maintains NLO accuracy for integrated (i.e. inclusive) quantities. 



Notice that the same formula ( |1.3D also describes the radiation of the hardest parton in MC@NLO, 
provided R" is identified with the shower Monte Carlo (i.e. HERWIG's) approximation of the real 
emission cross section. In fact, in MC@NLO two types of events are generated, called y and 
events. =5^ events correspond to the term proportional to B. In MC@NLO the corresponding un- 
derlying Bom kinematics is generated with a probability B{^B)d^B, while the hardest radiation 
kinematics is generated by the HERWIG shower algorithm. It was demonstrated in ref. ^ that the 
hardest radiation in HERWIG corresponds to the factor in square bracket multiplying B in eq. ([Qj). 



The J^f events correspond instead to the R—R'^ term in eq. ( |1.3| ). However, since R^ is now given by 
the HERWIG shower algorithm, there is no guarantee that the difference R — R^ should be positive, 
and this is why negative weighted events are an essential feature of MC@NLO. 

The fact that both in MC@NLO and POWHEG the hardest radiation can be described by a similar 
formula has allowed a better understanding of the agreement and discrepancies between the two 
approaches. First of all, one understands why most distributions compare very well in the two 
schemes (see for example ^ ^, 12]), since they are both described by a similar formula. A first 



area of discrepancy has emerged following the work of ref. [18]. In production, a dip in the 
rapidity distribution of the hardest jet of MC @NLO was found, that is not present neither in ALPGEN 
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Figure 1: Left: rapidity distribution of the hardest jet relative to the Higgs rapidity. The highest curve at 

y = is POWHEG+PYTHIA, next is POWHEG+HERWIG, then MCgNLO and then HERWIG alone. Right: 
MCSNLO, with B replaced by B in the generation of events. 



nor in POWHEG. It was shown later |l^ that this dip is a feature of MCSNLO that is present in 
several processes. It is particularly visible in Higgs production, as one can see from fig. The pure 
HERWIG result exhibits a dip at zero rapidity. MCgNLO has a similar dip, although relatively less 
severe. No dip is present in POWHEG. The origin of the dip is easily understood as a consequence of 
the presence of the dip in pure HERWIG and of formula ( |l.3| ). In fact, for large transverse momenta. 



formula ( |1.3| ) becomes equal to 

R'd<^. (1.8) 



do=Bx —d^ +[R- R'jd^ = Rd^ + 
B 



The second term is small in a perturbative sense, being of order a,. However, in the case of Higgs 
production, where NLO corrections amount to 100%, it is in fact of order 1. Since R" correspond to 
the HERWIG approximation to the real cross section, it has a dip that propagates into the MC@NLO 
cross section. We have checked this explicitly. If this argument is correct, by replacing the B cross 
section with B in MC @NLO, the dip should go away. This is in fact shown in fig. on the right plot. 

Another important discrepancy is found in the Higgs px distribution in the gg ^ h process, 
displayed in the left plot of fig. |2[ At high pr the POWHEG and MCgNLO spectra depart sensibly. 



The origin of this difference can again be understood starting from formula (1.8). In the POWHEG 
implementation of Higgs production, the choice R" = R was made, so that formula (1.8) yields 



B 

da = -xRd^. (1.9) 
B 

Because of the large Higgs production NLO corrections, the B/B factor, instead of being near 1 is 
close to 2. Thus NLO corrections amplify the whole pj spectrum in POWHEG. Of course, we could 
have chosen a different R" in order to reduce this effect, but in view of the plot on the right in fig. ^, 
this choice was maintained, since the px spectrum produced in this way is closer to the NNLO 
result. In MCgNLO this amplification works only for the y events, that become negligible at large 
Pt- Notice that in fig. ^ a mild kink is visible at 80-1 10 GeV in the MCgNLO curve, signalling the 
separation of the y and dominated regimes. 
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Figure 2: Left: rapidity distribution of the hardest jet relative to the Higgs rapidity. The highest curve at 
y = is POWHEG+PYTHIA, next is POWHEG+HERWIG, then MC@NLO and then HERWIG alone. Right: 
MCSNLO, with B replaced by B in the generation of y events. 



2. The POWHEG BOX 

In ref. the general formulation of POWHEG for complex processes was given, dealing with 
all complications having to do with the presence of several massless coloured lines, both in the 



final and in the initial state. In this case formula ( |1.3| ) should be applied to each radiating line. 
What a regular shower Monte Carlo does, starting from a Born process, is to let each coloured 
massless line radiate more light partons. Within POWHEG one does something similar. However, 
since in this case the radiation cross section is the exact cross section, rather than an Altarelli-Parisi 
approximation to it, we work the other way around, decomposing the full real cross section into 
a sum of terms, each of them singular in only one coUinear region. We then adopt a different 
phase space parametrization d^^sd^r for each singular region. After this is done, radiation 

is generated from each coUinear region, but only the hardest radiation is retained, and the event is 
generated according to its kinematics. This procedure reaches a considerable level of complexity, 
so we decided to implement it in a general way, in order to ease the POWHEG implementation of 



complex processes. This project, named the POWHEG BOX has been completed |2C], and will 



soon been published. It has been used to reproduce the previous implementation ofW/Z and 



single top production. Furthermore, it has been used to implement Z+jet production [13], and 



Higgs production in vector boson fusion []21[]. This last two processes are the first ones not already 
present in MC@NLO. 

Internally, the POWHEG BOX implements the NLO corrections using the FKS subtraction 



framework [g2|, g3|]. The author of an NLO calculation, willing to implement it into the POWHEG 
BOX framework, needs only to supply the following ingredients: a) the Born phase space; b) the 
hst of Born and Real processes in a specified format; c) the Born squared amplitudes B = |^^|, 
the colour correlated one B,y, and the spin correlated one Bj^.^r, for each Born subprocess (these 
ingredient are typically needed in NLO calculations performed with the subtraction method, and. 



for example, are generated automatically in the MadDipole package [24]); d) the real squared am- 
plitude, for all relevant partonic processes (this is easily obtained from programs like MadGraph 
[p5[]); e) the finite part of the virtual amplitude for all relevant partonic processes, which is the only 
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ingredient that may require dedicated work. The POWHEG BOX does all the rest, including the 
automatic implementation of the FKS subtractions. 



3. Conclusions 



Much progress has been achieved in POWHEG in this past year. Several features of the method 
have been better understood in simple terms, and its comparison with MC@NLO has been 
considerably clarified. An important step towards automation has been achieved with the 
construction of the POWHEG BOX. Thanks to this framework, a relatively complex process like 
VBF Higgs production has been implemented in a very short time. It is our hope that the POWHEG 
BOX will allow authors of NLO calculations to implement their work in the POWHEG framework 
with a limited effort. 
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